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^ ; ABSTRACT 

I ^ . We study the gamma-ray emissions from an outer-magnetospheric potential gap 

^ ! around a rotating neutron star. Migratory electrons and positrons are accelerated by 

CSJ I the electric field in the gap to radiate copious gamma-rays via curvature process. Some 

^ of these gamma-rays materialize as pairs by colliding with the X-rays in the gap, leading 

^ . to a pair production cascade. Imposing the closure condition that a single pair produces 
one pair in the gap on average, we explicitly solve the strength of the acceleration field 

^ I and demonstrate how the peak energy and the luminosity of the curvature-radiated, 

l/^ ■ GeV photons depend on the strength of the surface blackbody and the power-law emis- 

O , sions. Some predictions on the GeV emission from twelve rotation-powered pulsars are 

Q i presented. We further demonstrate that the expected pulsed TeV fiuxes are consistent 

^ I with their observational upper limits. An implication of high-energy pulse phase width 

P^' versus pulsar age, spin, and magnetic moment is discussed. 

O 
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1. Introduction 

The EGRET experiment on the Compton Gamma 
Ray Observatory has detected pulsed signals from 
seven rotation-powered pulsars (e.g., Nolan et al. 
1996, and references therein; Kapsi et al. 2000): 
Crab, Vela, Geminga, PSRB1706-44, PSRB1951+32, 
PSR B1046-58, and PSR B1055-52, with PSR B0656+14 
being a possible detection (Ramanamurthy et al. 
1996). The modulation of the 7-ray light curves at 
GeV energies testifies to the production of 7-ray radi- 
ation in the pulsar magnetospheres either at the polar 
cap (Harding, Tademaru, & Esposito 1978; Daugh- 
erty & Harding 1982, 1996; Dermer & Sturner 1994; 
Sturner, Dermer, & Michel 1995; Shibata, Miyazaki, 
& Takahara 1998; Miyazaki & Takahara 1997; also 
see Scharlemann, Arons, & Fawley 1978 for the slot 
gap model), or at the vacuum gaps in the outer mag- 
netosphere (Chen, Ho, & Ruderman 1986a, b, here- 
after CHR; Chiang & Romani 1992, 1994; Romani 
and Yadigaroglu 1995; Romani 1996; Zhang & Cheng 
1998; Cheng & Zhang 1999). Effective 7-ray produc- 
tion in a pulsar magnetosphere may be extended to 
the very high energy (VHE) region above 100 GeV as 
well; however, the predictions of fluxes by the current 
models of 7-ray pulsars are not sufficiently conclu- 
sive (e.g., Cheng 1994). Whether or not the spectra 
of 7-ray pulsars continue up to the VHE region is a 
question which remains one of the interesting issues 
of high-energy astrophysics. 

In the VHE region, positive detections of radiation 
at a high confidence level have been reported from 
the direction of the Crab, B1706-44, and Vela pulsars 
(Bowden et al. 1993; Nel et al. 1993; Edwards et 
al. 1994; Yoshikoshi et al. 1997; see also Kifune 1996 
for a review), by virtue of the technique of imaging 
Cerenkov light from extensive air showers. However, 
with respect to pulsed TeV emissions, only the upper 
limits have been, as a rule, obtained from these pul- 
sars (see the references cited just above). If the VHE 
emission originates the pulsar magnetosphere, rather 
than the extended nebula, significant fraction of them 
can be expected to show a pulsation. Therefore, the 
lack of pulsed TeV emissions provides a severe con- 
straint on the modeling of particle acceleration zones 
in a pulsar magnetosphere. 

In fact, in CHR picture, the magnetosphere should 
be optically thick for pair production in order to re- 
duce the TeV flux to an unobserved level by absorp- 
tion. This in turn requires very high luminosities of 



tertiary photons in the infrared energy range. How- 
ever, the required fluxes are generally orders of mag- 
nitude larger than the observed values (Usov 1994). 
We are therefore motivated by the need to contrive 
an outer gap model which produces less TeV emis- 
sion with a moderate infrared luminosity. 

High-energy emission from a pulsar magnetosphere, 
in fact, crucially depends on the acceleration electric 
fleld, along the magnetic fleld lines. It was Hi- 
rotani & Shibata (1999a, b; hereafter Paper I, II) who 
flrst solved the spatial distribution of iSy together 
with particle and 7-ray distribution functions. They 
explicitly demonstrated that there is a stationary so- 
lution for an outer gap which is formed around the 
null surface at which the local Goldreich- Julian charge 
density 

27rc[l - (nw/c)^] ^ ' 

vanishes, where is the component of the magnetic 
fleld along the rotation axis, fl refers to the angular 
frequency of the neutron star, m indicates the dis- 
tance of the point from the rotation axis, and c is 
the speed of light. Subsequently, Hirotani (2000a, 
hereafter Paper IV) investigated the 7-ray emission 
from an outer gap, by imposing a gap closure condi- 
tion that a single pair produces one pair in the gap 
on average. He demonstrated that E^^ becomes typi- 
cally less than 10% of the value assumed in CHR and 
that the resultant TeV flux is sufficiently less than 
the observational upper limit of the pulsed fiux, if the 
outer gap is immersed in a X-ray field supplied by the 
blackbody radiation from the whole neutron star sur- 
face and/or from the heated polar caps. In this paper, 
we develop his method to the case when a magneto- 
spheric power-law component contributes in addition 
to the blackbody components. 

In the next section, we formulate the gap closure 
condition. Solving the condition in §3, we investigate 
the acceleration field and the resultant 7-ray emis- 
sions as a function of the X-ray field. In §4, we further 
apply the theory to twelve rotation-powered pulsars 
and predict the absolute fiuxes of TeV emission from 
their outer gaps. In the final section, we discuss the 
validity of assumptions and give some implications on 
pulse profiles of GeV emissions. 

2. Electrodynamic structure of the gap 

We first describe the acceleration field in §2.1, then 
consider the energy of curvature-radiated 7-rays in 
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§2.2, the X-ray field illuminating the gap in §2.3 and 
the pair production mean free path in §2.4. We fur- 
ther formulate in §2.5 the gap closure condition that 
one of the copious 7-rays emitted by a single pair ma- 
terialize as a pair in the gap on average. We finally 
present the resultant 7-ray properties in §§2.6 and 2.7. 

2.1. Acceleration Field in the Gap 

In this paper, we consider an outer gap in the 
following rectilinear coordinate: x is an outwardly 
increasing coordinate along the magnetic field lines, 
while z is parallel to the rotational axis. We define 
a; = to be the intersection between the last open 
field line and the null surface where 3^ = (fig. |^). 
If we assume that the transfield thickness of the gap 
is larger than the gap width along the field lines, we 
can neglect dz'^ term compared with dx^ one in the 
Poisson equation for the non-corotational potential, 
Then we can Taylor-expand pcj around a; = to 
obtain 

- ^ = -4^Ax, (2) 

where A is the expansion coefficient of pcj at x = 
0. Since the toroidal current flowing near the light 
cylinder is unknown, we simply approximate with 
its Newtonian value. Then A is given by 

~ 271^0"* 1 - (Jlro sin 6*0/0)2 
■ sin 26*0 cos(6'o — ai) -I- cos 26*0 sin(6'o — ai) 



(3) 



where refers to the magnetic dipole moment of the 
neutron star; tq is the distance of the gap center (a; = 
0) from the neutron star, and Bq is the polar angle 
of the gap center (in the first quadrant). They are 
related with a\ as 

sin2(6'o — ai) 



w-LQ sin^LC sin2(0LC — ai) 



(4) 



and 



3 tan Q!i + \/9 tan^ ai + 8 
tant^o = -z , (5) 



where the light-cylinder radius is defined by 

G7LC = ^ = 3.0 X lO^^Oa^^cm (6) 

and the colatitude angle 0lc at which the last-open 
fieldline intersects with the light cylinder is implicitly 
solved from 



^2 = rJ/lO^rad % ^ . In the order of magnitude, 
A ^ Pg.t/'^lc holds. Its exact values are 3.19 x 
lO^i^f^s V30 esu for a.-, = 30° and 1.22 x 10-"f72 Vso 
esu for Qfj = 45°, where /iao = /i/10'^"G cm'^. 

Integrating equation (|^), wc obtain the accelera- 
tion field E\\(x) = -d<^/dx = E\\{0) - 2ttAx^, where 
i?l|(0) refers to the value of at a; = 0. Defining 
the boundaries of the gap to be the places where iSy 
vanishes, we obtain i?||(0) = 2Tr AH^. It is notewor- 
thy that the non- vanishment of dE\\/dx at a; = ±H is 
consistent with the stability condition at the plasma- 
vacuum interface if the electrically supported magne- 
tospheric plasma is completely-charge-separated, i.e., 
if the plasma cloud &i x < —H \s, composed of elec- 
trons alone (Krause-Polstorff & Michel 1985a, b). We 
assume that the Goldrcich- Julian plasma gap bound- 
ary is stable with = on the boundary, x — ^H. 

We can now evaluate the typical strength of by 
averaging its values throughout the gap as follows: 



En 



H 



1 

H 



dx {En{0) - 2ttAx'^) 



(8) 



We shall use this E\\ as a representative value of 
the acceleration field in the gap. Defining the non- 
dimensional gap width h = H/w-^c and substituting 
the values of A, we obtain 



3 2 V 



(9) 



where Ce = 3.61 x lO^o for a-, = 30°, and 1.39 x 10" 
for ai = 45°. 

The voltage drop in the gap is given by 



Fgap = Ell ■ = s'"^^ 



(10) 



2 cot(6'LC - ai) sin 6'lc + cos 6'lc = 0; 



(7) 



For an outer gap which extends to the light cylinder 
{H ~ 0.5n7Lc), V^ap becomes as large as the available 
electromotive force (EMF) exerted on the spinning 
neutron star surface, K ~ 10^^(02/0.5)2/i3oV. 

2.2. Energy of curvature-radiated 7-rays 

The most effective assumption for the particle mo- 
tion in the gap arises from the fact that the velocity 
saturates immediately after their birth in the balance 
between the radiation reaction force and the electric 
force. The reaction force is mainly due to curvature 
radiation if the gap is immersed in a moderate X-ray 
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field. Equating the electric force eE\\ and the radi- 
ation reaction force, we obtain the Lorentz factor of 
saturated particles as follows: 



r = 



2e 



1/4 



(11) 



where i?c and e refers to the curvature radius of the 
magnetic field line and the magnitude of the charge 
on the electron. T is related with h = H/wi^c by 



(12) 



where Cr = 9.59 x lO'^ for a-, = 30°, and 1.34 x 10^ 
for ai = 45°; R0.5 = i?c/(0.5n7Lc). 

However, if the gap width H is less than the length 
scale for the particles to be accelerated up to F, the 
typical Lorentz factor should be rather estimated by 
eE^^H/nicC^, because typical particles are accelerated 
by the potential V^ap/2 = Therefore, we can 

evaluate the Lorentz factor as 



r = niin(r, e£^||-ff/moC^) 



(13) 



Using this F, we obtain the central energy of cur- 
vature radiation. 



3F3 he 



Er = — — 



(14) 



where h is the Planck constant divided by 27r. In 
this paper, we adopt the gray approximation in the 
sense that all the 7-rays are radiated at energy E^. In 
the final section of Paper IV, we demonstrated that 
this gray approximation gave a good estimate of the 
gap width and other quantities describing the gap, 
by comparing them with those obtained in the non- 
gray cases in which the Boltzmann equations of parti- 
cles and 7-rays were solved together with the Poisson 
equation for <i>. 

2.3. X-ray field 

Before proceeding to the discussion on pair pro- 
duction mean free path, it is desirable to describe the 
X-ray field that illuminates the outer gap. X-ray field 
of a rotation-powered neutron star within the light 
cylinder can be attributed to the following three emis- 
sion processes: 

(1) Photospheric emission from the whole surface of a 
cooling neutron star (Greenstein & Hartke 1983; Ro- 
mani 1987; Shibanov et al. 1992; Pavlov et al. 1994; 



Zavlin et al. 1995). 

(2) Thermal emission from the neutron star's polar 
caps which are heated by the bombardment of rela- 
tivistic particles streaming back to the surface from 
the magnetosphere (Kundt & Schaaf 1993; Zavlin, 
Shibanov, & Pavlov 1995; Gil & Krawczyk 1996). 

(3) Non-thermal emission from relativistic particles 
accelerated in the pulsar magnetosphere (Ochelkov 
& Usov 1980a,b; El-Gowhari & Arponen 1972; As- 
chenbach & Brinkmann 1974; Hardee & Rose 1974; 
Daishido 1975). 

The spectrum of the first component is expected 
to be expressed with a modified blackbody. However, 
for simplicity, we approximate it in terms of a Plank 
function with temperature fcTg, because the X-ray 
spectrum is occasionally fitted by a simple blackbody 
spectrum. We regard a blackbody component as the 
first one if its observed emitting area is comparable 
with the whole surface of a neutron star, A* = Anr^:'^, 
where denotes the neutron star radius. We take ac- 
count of both the pulsed and the non-pulsed surface 
blackbody emission as this soft blackbody component. 

As for the second component, we regard a black- 
body component as the heated polar cap emission 
if its observed emitting area is much smaller than 
A^. We approximate its spectrum by a Planck func- 
tion. We take account of both the pulsed and the 
non-pulsed polar cap emission as this hard blackbody 
component. 

Unlike the first and the second components, a 
power-law component is usually dominated by a neb- 
ula emission. To get rid of the nebula emission, which 
illuminates the outer gap inefficiently, we adopt only 
the pulsed components of a power-law emission as the 
third component. 

2.4. Pair Production Mean Free Path 

In this subsection, we draw attention to how the 
pair production mean free path is related with the 
X-ray field described in the previous section. To this 
end, we first give the threshold energy for X-ray pho- 
tons to materialize. Then, we consider the mean free 
path for a 7-ray photon to materialize as a pair in 
a collision with one of the soft blackbody X-rays in 



2.4.2, the hard blackbody ones in § 2.4.3, and the 



magnetospheric, power-law X-rays in § 2.4.4. We fi- 



nally summarize how the real mean free path should 
be computed under the existence of these three X-ray 
components in § 2.4.5. 
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2.4-1- Threshold Energy 

The threshold energy for a soft photon to mate- 
riaUze as a pair in a coUision with the 7-ray having 
energy meC^e^, = E^, is given by 



E. 



th 



1 — cos 0ab cos ( 



(15) 



where cos 0ab cos 9c refers to the cosine of the three- 
dimensional coUisional angle between the X-ray and 
the 7-ray. The lower bould of 2/(1 — cos0abCOS0c) 
is unity, which is realized if the two photons head-on 
collide. 

To evaluate cos (^ab cos 0c, we must consider 7- 
ray's toroidal momenta due to aberration. At the 
gap center, the aberration angle 0ab is given by 
sin~^(ro sin6'o/'Ci7Lc). In the case of ai — 30°, we 
obtain ^ab = 20.4°. The collisional angle on the 
poloidal plane becomes 9c = 90° —6*0 = 21.6° (or 6'c = 
90° + 9o = 158.4°) for outwardly (or inwardly) prop- 
agating 7-rays, where is computed from equation 
(|). Therefore, we obtain cos^abcos^c = ±0.855, 
where the upper and the lower sign correspond to the 
outwardly and inwardly propagating 7-rays, respec- 
tively. On the other hand, in the case of — 45°, we 
have 0ab = 14.0° and 6'c = 15.6° (or 9c = 164.3°) to 
obtain cos</>ab cos 6'c = ±0.923. 

2-4-2- Soft Blackbody Component 

We first consider the pair production mean free 
path. As, for a 7-ray photon to materialize in a col- 
lision with one of the soft blackbody X-rays. In a 
realistic outer gap, both the outwardly and inwardly 
propagating 7-rays contribute for Ag. Therefore, As 
is given by an arithmetic average of these two con- 
tributions as follows (Paper IV; see also eq. [3.1] in 
Blandford & Levinson 1995 for the factor 1 ± /Zc): 



1 

aT 



re(l - ^c) 



a€x— — CTp(e^,ex,/ic) 



f°° dNs 
±(1 - k)(1 + /ic) / de^- — f7p(e^,ex,-/.tc), 

(16) 

where fic = | cos(/)ab cos^d; the pair production cross 
section is given by (Berestetskii et al. 1989) 



o'p(e^,ex,/ic) 
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crT(l — v'^) 



(3- w^)ln 



1-t) 



- 2v{2-v^ 



1 



(17) 



y 1 - Mc <^t<^^ 

where ctt is the Thomson cross section, and Cx = 
Ex/rUcC^ refers to the dimensionless energy of an 
X-ray photon. We may notice here that the non- 
dimensional threshold energy {Eth/mcC^) appears in 
the lower bound of the integral in equation (|l6|). 

When a blackbody emission dominates the X-ray 
field, H occasionally becomes as large as 0.3ci7LC- 
In this case, we must take the dilution effect of the 
surface radiation into account. The number den- 
sity of the soft X-rays between energies TOcC^Cx and 
moC^(ex + rfcx) at position x is given by 



dN, 



1 + 2sin^ 



UN, 

\ *x 



(18) 

Here, the number density at the gap center is given 
by the Planck law 

3 



dN, 



1 

4^ 



m,cC 
cfi 



f A, 



V47rr^y exp(ex/As) - 1' 
(19) 

where A, indicates the observed emission area of the 
soft blackbody; Ag is defined by 

As^-^, (20) 
mcC 

where kTg refers to the soft blackbody temperature 
measured by a distant observer. Since the outer gap 
is located outside of the deep gravitational potential 
well of the neutron star, the photon energy there is 
essentially the same as the distant observer measures. 

The first (or the second) term in equation ( p^ rep- 
resents the contribution from the outwardly (or in- 
wardly) propagating 7-rays; For ai = 30°, we adopt 
fic = 0.855, whereas for 45°, /ic = 0.923. The weight 
K reflects the ratio of the fluxes between the outwardly 
and inwardly propagating 7-rays. For example, if k 
were to be 1.0, only outwardly propagating 7-rays 
would contribute for the pair production. From Pa- 
per III, we know that the flux of the outwardly propa- 
gating 7-rays is typically about ten times larger than 
the inwardly propagating ones. In what follows, we 
thus adopt K = 0.9. This value of k was justified in 
the final section of Paper IV. 

2-4-3- Hard Blackbody Component 

In this subsection, let us consider the case when 
the X-ray field is dominated by the hard blackbody 
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component. In the same manner as the soft blackbody 
component, the hard blackbody component gives the 
following mean free path for pair production: 



-^h J2/(1-Mc)e-, 
+ (l-K)(l+/ic) 



acx-i — CTp^e-y, Ex, lie) 

CLEy- 



2/(l+Mo)£ 



wEx-i — crp^e^, ex, — A^cj, 
(21) 



where //c takes the same value as the soft-blackbody- 
dominated case and k = 0.9. The number density of 
the hard blackbody X-rays between energies me(?e-i^ 
and moC^(ex + dcx) at position x is given by 



dfx 

f dN^ 
Wcx / 



1 + 2 sin 6*0- 



X 



( ^ 



1 / mr C 
4^ 



,2x 3 



Ah 



(22) 



47rr5y exp(ex/Ah) - 1 
(23) 

where Ah denotes the observed emission area of the 
hard-blackbody emission. Ah is defined by 



Ah 



(24) 



where fcTh refers to the hard blackbody temperature 
measured by a distant observer. 

2.4-4- Power-law Component 

Since the secondary photons emitted outside of the 
gap via synchrotron process will be beamed in the 
same direction of the primary 7-rays, the typical col- 
lision angle on the poloidal plane can be approximated 
ash = H/wlc rad. It follows that the mean free path 
corresponding to the power-law emission becomes 



1 



where 



^ = (l-^c)/ dex^^^crp(e^,ex,Afc), (25) 



2 



1 - 



(26) 



dNp\/dex refers to the number density of the power- 
law X-rays between energies rUcC^e^ and mcC^{e^ -\- 
de~f); we adopt 



H 



To sm ^0 



fab 



cos h 



(27) 



for the power-law component. We may notice here 
that both the 7-rays and the magnetospheric, power- 
law X-rays suffer aberration and that the resultanlt 
collision angle in the azimuthal direction is less than 
4>ab and can be assumed to be (i?/ro sin 6'o)0ab- The 
dependence of dN-pi/de^ on x is unknown. We thus 
simply assume that it is constant throughout the gap 
and specify the form as 



dN, 



pi 



dey, 



= Afpiex" (e. 



< Cx < En 



,). (28) 



The photon index a is typically between —2 and — 1 
for a pulsed, power-law X-ray component in hard X- 
ray band (e.g., Saito 1998). 

2.4.5. Pair Production Mean Free Path 

The true mean path, Ap, is determined by the 
component that dominates the X-ray field, or equiva- 
lently, the smallest one among As, Ah, and Api. There- 
fore, we can reasonably evaluate Ap as 



1 

a: 



1 

aT 



1 



1 

Api 



(29) 



When a pulsar is young, the third term will dom- 
inate because of its strong magnetospheric emission. 
As the pulsar evolves, the first term will become dom- 
inant owing to the diminishing magnetospheric emis- 
sion. As the pulsar evolves further, polar cap heat- 
ing due to particle bombardment begins to dominate; 
therefore, the second term becomes important. 

2.5. Gap Closure 

The gap width 2H is adjusted so that a single pair 
produces copious 7-ray photons (of number iV^) one 
of which materializes as a pair on average. Since a 
typical 7-ray photon runs the length H in the gap 
before escaping from either of the boundaries, the 
probability of a 7-ray photon to materialize within 
the gap, Nj^^, must coincide with the optical depth 
for absorption, H/Xp. Considering the position de- 
pendence of Ap on X, we obtain 



where 



1 



1 dx 

2 J-H '^P 



(30) 



(31) 



Here, Ap is given by equation (g9|). Equation (|3l|) is 
derived as follows: A single e"*" or e~ emits 7-rays at 
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the rate P = 2ce^V'^ / {ZR^^) . Dividing P with the 
averaged photon energy, 'iT^hc/{2Rc), we can esti- 
mate the number of 7-rays emitted per unit time by 
a single e+ or e~, N^. Noting that a typical parti- 
cle runs length H in the gap on average, we obtain 



= {H/c)Nj, which reduces to equation (31) 



Substituting equations (nSl), and into 



(30), we obtain the gap closure condition 



1 



H 

As,0 



H 



h,0 



I{h) 



H 

Xpi ' 



(32) 



where As,o and Ah,o refer to the values of As and Ah 
at a; = 0, and 



lih) 



tan ^ (tan 9o + h sec ^o) — tan ^ (tan do ~ h sec 6*0) 



2h cos 00 

(33) 

For a very thin gap {h ~ H/vji^q ^ 1), I{h) ap- 
proaches to unity. When the surface blackbody com- 
ponents dominate the X-ray field, h can be as large as 
0.3; as a result, the r^^ variation of the X-ray density 
becomes important. We take account of this effect in 
function I{h). When the power-law component dom- 
inates, on the other hand, /i <C 1 is usually satisfied 
(see discussion in § O); therefore, we only evaluate 
the X-ray density at the gap center in equation (|3^). 
Combining equations (31) and (^2|), and represent- 
ing r, — Ec/mcC^ in terms oi h = _ff/n7Lc, we 
finally obtain the equation describing ft, as a function 
of B5, il2, kTs, As, kT\i, and Ah. Once h is solved, 
other quantities such as , F, Ec can be computed 
straightforwardly. 

2.6. Luminosity of GeV emissions 

Let us first consider the luminosities of curvature- 
radiated 7-rays. The luminosity, Lq^y, can be esti- 
mated by multiplying the total number of positrons 
and electrons in the gap, N^, the number of 7-rays 
emitted per particle per unit time, / (H / c) , and 7- 
ray energy, Ec . Supposing the conserved current den- 
sity is xcpGJ = x(^-B/27re), and assuming that the 
gap exists tt{H/tui^c) rad in the azimuthal direction, 
we obtain 

~ -^-ro sm0o ■ ■ 2H. (34) 

/Tree wlc 

where D± is the transfield thickness of the gap on the 
poloidal plane. The distance of the center of the gap 
from the rotation axis, rpsin^Oi becomes 0.393ti7Lc 
for ai = 30°, and 0.285n7Lc for a-, = 45°. The 



strength of magnetic field at the gap center can be 
given by 



B = 3.70 X lOVl + 3cos2(6'o - ai) 



(35) 



As discussed in §2.4 in Paper IV, we adopt x — 0-1 
and D± — 0.3ci7LC as typical values in this paper. 
From equations ( pT] ) and ([T^), we obtain 

2/ \ el7Brosin6lo 4 2 .orn 

^Gcv = - 1 X— — ) ,n — h , (36) 



3 V ^LC/ (i?c/wLc)^ 



where F is given by equation ([1^) . If F is saturated, 
equation (^ or (|l2|) gives 

LgcV ^ CgcV f x— ) f^2'*M30^/i'*ergs s~\ (37) 

where CgcV = 1.19x lO^^ for a; = 30°, and 1.20x10'*° 
for ai = 45°. 

2.7. Luminosity of TeV emissions 

The relativistic particles produce 7-rays mainly via 
curvature radiation as described in the preceding sec- 
tions. However, even though energetically negligible, 
it is useful to draw attention to the TeV 7-rays pro- 
duced via inverse-Compton (IC) scatterings. Since 
the particles' Lorentz factor becomes ~ 10^'^, it is 
the infrared photons with energy 0.01 eV that con- 
tribute most effectively as the target photons of IC 
scatterings. Neither the higher energy photons like 
surface blackbody X-rays nor the lower energy pho- 
tons like polar-cap radio emission contribute as the 
target photons, because they have either too small 
cross sections or too small energy transfer when they 
are scattered. On these grounds, we obtain the follow- 
ing upper bound for the luminosity of the IC-scattered 
7-rays: 



XtcV < 5.63 X lO^^BgLgo ( sin 6*0 

F/i^ergs s~ 



X — 
(38) 



where L30 refers to the luminosity of infrared photons 
that can be scattered up to TeV energy range in the 
unit of lO^^ergs s'^ B5 = B/{W^G). The inequality 
comes from the fact that the scattered 7-rays cannot 
have energies greater than TrricC^. 
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The vF^ flux of IC-scattered 7-rays, {vF,j)tcV 
[Jy Hz], can be readily computed as 



TcV 



10 



-20 



L 



ToV 



era's s 



1 ster 



di~'Jy Hz, 
(39) 

where wxeV refers to the sohd angle in which the TeV 
7-rays are emitted, and di = d/lkpc. 

An analogous relation holds for an infrared flux, 
(i/i^y)o.oieV: if the infrared luminosity, -Lq.oioV, is 
emitted in a solid angle wq.oioV ster. We thus obtain 
the flux ratio 



ToV 



5.63 X 10 



-5 f ^LC . 



sin 9q 



O.OloV 



xrlx^^TT^^^^^/i^. (40) 



V '^0 

iulc ) lO^G wtcV 

Since we take the ratio of vF^, flux at two different 
energies, the uncertainties arising from the distance 
disappears on the right-hand side. It follows from 
equation ( ^o|) that the ratio becomes at most 10'^ in 
order of magnitude, because x(Z3j^/ti7Lc) < 1, /i < 
0.5, and T ^ 10^'^ hold in general. 

3. X-ray field vs. 7-ray emission 

Before proceeding to an application to individual 
pulsars, it will be useful to investigate the general 
properties of 7-ray emission as a function of the X- 
ray fleld illuminating the gap. To this aim, we first 
demonstrate how the gap width h = H/wi,c depends 
on the X-ray field in §3.1, by solving the gap clo- 
sure condition ( |3^ . We then present the energies and 
the luminosities of the curvature-radiated and the IC- 
scattered 7-rays in §3.2. Throughout this section, we 
adopt fJ2 = 0.5, /i3o = 3.0, = = AirrJ^, kT^ = 
200eV, a = -1.5, e^in = 0.1keV/511keV = 0.00020, 
and e^ax = 100keV/511keV 0.20. 

3.1. Gap width 

The results of the gap half width divided by the 
light cylinder radius are presented in figure ^j. The ab- 
scissa is kTs in eV. The model parameters of the X-ray 
field for the six curves are summarized in table 1. For 
the three thick curves, the hard blackbody component 
is not included (i.e., A\^ = Q). The thick solid, dashed, 
and dotted curves correspond to A^pi = 0, 10^", 10^^, 
respectively. Therefore, the thick solid curve corre- 
sponds to the least dense X-ray field thereby gives 
the greatest h = H/wi,c for a specific value of kTs- 



For the three thin curves, on the other hand, A^pi is 
set to be 0, so that the power-law component does 
not contribute. The thin solid, dashed, and dotted 
curves correspond to AhM* = 10^"*, lO^'^, 10~^, re- 
spectively. 

First of all, it follows from the figure that H (or 
equivalently h for a fixed fi) is a decreasing function 
of kTs- The reason is as follows: If kTs increases, the 
number density of target soft photons above thresh- 
old energy Ns{E > Eth) increases for a fixed value of 
Eth- The increased Ns{E > Eth) results in the de- 
crease of Ap , which reduces H (eq. Q ) . Accurately 
speaking, the reduced H results in a decrease of 
and hence E^, thereby increases Eth- The increased 
Eth decreases Ns{E > Eth) to partially cancel the 
initial decrease of Ap. In addition, the reduction of 
H implies the reduction of the emitting length for a 
particle, thereby decreases Nj and partially cancel the 
initial decrease of H (see eq. [^). Nevertheless, both 
of the two effects are passive; therefore, the nature of 
the decrease of H with increasing kTs is unchanged. 

The second thing to note is that h decreases with 
increasing ^h, as indicated by the three thin curves 
in figure |[ This is because when Ah increases the 
number density of the hard blackbody component, 
Nh{Ex > Eth), increases. This in turn leads to the 
decrease of Ap, which results in the decrease of h. If 
An is as small as 10^*A^, (the thin solid line), the X- 
ray field is dominated by the power-law component 
only in the small kTs range; this can be understood 
because the thin solid line significantly deviates from 
the thick solid line at kTs < 30eV. However, if Ah is 
as large as 10~^ (the thin dotted line), the X-ray field 
is dominated by the hard blackbody component up to 
very high kT^ (- lOOeV). 

The third thing is that h decreases with increasing 
A^pi. For a very strong magnetospheric emission (the 
thick dotted line), h becomes not more than 0.1. 

In short, the gap width is a decreasing function of 
of the X-ray number density, regardless of the com- 
ponent that dominates the X-ray field. 

3.2. Gamma-ray Luminosities 

Let us now consider the luminosities of the curvature- 
radiated 7-rays. Substituting the results of h into 
(^), we obtain Lgcv; in figure ^, we present the re- 
sults for X = 0.1 and D± = O.Stulc- It follows from 
this figure that XgcV is a decreasing function of the 
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Table 1: Model parameters of X-ray field 









v4hM* 




—a 




-^max 






cV 




cm-'^ 




keV 


keV 


thick solid 


1.0 


200 








2.0 


0.1 


100 


thick dashed 


1.0 


200 





IQio 


2.0 


0.1 


100 


thick dotted 


1.0 


200 





IQie 


2.0 


0.1 


100 


thin solid 


1.0 


200 


10-4 





2.0 


0.1 


100 


thin dashed 


1.0 


200 


10-3 





2.0 


0.1 


100 


thin dotted 


1.0 


200 


10-2 





2.0 


0.1 


100 



X-ray number density, regardless of the component 
that dominates the X-ray field. This is because the 
increase of the target X-ray photons results in the 
decreases of H. 

We next present the expected ratio of the TeV 
and the infrared fluxes in figure U, by using equation 
(40). It follows that the i^Fi, flux of TeV emission 
will not exceed lO^^Jy Hz for a typical infrared flux 
(< 10^ Jy Hz). Therefore, we can conclude that the 
pulsed TeV emission from rotation-powered pulsars 
cannot be detected in general by the current ground- 
based telescopes. 

4. Application to Individual Pulsars 

In this section, we apply the theory to the eleven 
rotation-powered pulsars of which X-ray field at the 
outer gap can be deduced from observations. We first 
describe their X-ray and infrared fields in the next 
two subsections, and present resultant GeV and TeV 
emissions from individual pulsars in § 4.3 and 4.4. 

4.1. Input X-ray Field 

We present the observed X-ray properties of indi- 
vidual pulsars in order of spin-down luminosity, iSrot 
(table 2). We assume emin = 0.1keV/511keV and 
Emax = 100keV/511keV for homogeneous discussion. 
Crab From HEAO 1 observations, its X-ray field 
is expressed by a power-law with a = — 1.81 ± 0.02 
in the primary pulse (PI) phase (Knight 1982). From 
the nebula and background subtracted counting rates, 
its normalization factor of this power-law emission be- 
comes A^pi = 5.3 X 10^^(f{ro/wLc)^^, where d refers 
to the distance in kpc. 

B0540-69 From ASCA observations in 2-10 keV 
band, its X-ray radiation is known to be well fitted 
by a power-law with a = —2.0. The unabsorbed lumi- 
nosity in this energy range is 1.3 x lO'^^erg s""'^, which 



leads to A^pi = 9.0 x I0^d^{ro/wi,c)~^ (Saito 1998). 
B1509-58 From ASCA observations in 2-10 keV 
band, its pulsed emission can be fitted by a power- 
law with a = —1.1. The unabsorbed flux in this en- 
ergy range is 2.9 x lO-^^erg/cm^/s, which leads to 
iVpi = 9.3 X 10"d2(ro/n7Lc)"'. 

J1617-5055 and J0822-4300 These two pulsars 
have resemble parameters such as f2 ~ QOrads"^, 
/i ^ 10'^*''®G cm^, and characteristic age r ~ 8 x 
lO^yr. From the ASCA observations of J1617-5055 
in 3.5-10 keV band, its pulsed emission subtracted 
the background and the steady components can be 
fitted by a power-law with a — —1.6 (Torii et al. 
1998). Adopting the distance to be 3.3 kpc (Coswell 
et al. 1975), we can calculate its unabsorbed flux 
as 3.1 X lO'i^erg/cm^/s, which yields A^pi = 6.5 x 
10i*'(i2(ro/n7Lc)^^- On the other hand, the distance 
of J0822-4300 was estimated from VLA observations 
as d = 2.2±0.3kpc (Reynoso et al. 1995). ROSAT ob- 
servations revealed that the soft X-ray emission of this 
pulsar is consistent with a single-temperature black- 
body model with kT^ = 0.28 ± O.lOkeV and ~ 
0.04A,(d/2.2)2 (Petre, Becker, & Winkler 1996). 
Vela From ROSAT observations in 0.06-2.4 keV, 
the spectrum of its point source (presumably the pul- 
sar) emission is expressed by two components: Sur- 
face blackbody component with fcTg — 150eV and 
As = 0.066A,(d/0.5)2, and a power-law component 
with a — —3.3. However, the latter component does 
not show pulsations; therefore, we consider only the 
former component as the X-ray field illuminating the 
outer gap. 

B1951+32 From ROSAT observations in 0.1-2.4 
keV, the spectrum of its point source (presumably the 
pulsar) emission can be fitted by a single power-law 
component with a — —1.6 and intrinsic luminosity 
of 2.3 X 1033(d/2.5)2erg s-^ (Safi-Harb & Ogelman 
1995), which yields iVpi = 9.1 x 10"d2(ro/n7Lc)"^- 
The extension of this power-law is consistent with the 
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Table 2: Input X-ray field 



pulsar 


distance 




logio A* 















a 




kpc 


rad s 


lg(G cm^) 


eV 




eV 




cm 






Crab 


2.49 


188.1 


30.53 












1 


8 


B0540-69 


49.4 


124.7 


31.00 










1014.15 


2 





B1509-58 


4.40 


41.7 


31.19 










1014.04 


1 


1 


J1617-5055 


3.30 


90.6 


30.78 










1012.64 


1 


6 


J0822-4300 


2.20 


83.4 


30.53 


280 


0.040 












Vela 


0.50 


61.3 


30.53 


150 


0.066 












B1951+32 


2.5 


159 


29.68 










1013.55 


1 


6 


B1821-24 


5.1 


2060 


27.35 










1016.36 


1 


2 


B0656+14 


0.76 


15.3 


30.67 


67 


4.5 


129 


10-1.49 


105.25 


1 


5 


Geminga 


0.16 


26.5 


30.21 


48 


0.16 






105.00 


1 


6 


B1055-52 


1.53 


31.9 


30.03 


68 


7.3 


320 


10-3.64 








J0437-4715 


0.180 


1092 


26.50 


22 


0.16 


95 


10-3.28 










Fig. 2. — Examples of the gap half width divided by 
light cylinder radius as a function of kTg [eV]. The 
Fig. 1.— A side view of a hypothetical outer magne- j^put parameters of the X-ray field of each curve are 

tospheric gap in which pair production cascade takes i^g^^^ table 1. For all the six curves, pulsar param- 

Place. eters are fixed as f2 = 50 rad s'^ fi = W^°G cm^ 

and ai = 30°. 
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kT [eV] 

Fig. 3. — Luminosity [ergs s~^] of the curvature- 
radiated 7-rays. The abscissa represents fcTg in cV 
unit. The six curves correspond to the same parame- 
ter set as in figure 2 (or table 1). 




kT [eV] 

Fig. 4. — Ratio of the fluxes between TeV and in- 
frared energy ranges. The abscissa represents kTs in 
eV unit. The TeV 7-rays are radiated by the primary 
positrons and electrons in the gap via IC scatterings, 
while the infrared photons are presumably emitted 
by tertiary pairs outside of the gap via synchrotron 
process. The six curves correspond to the same pa- 
rameter set as in figure 2 (or table 1). 



upper limit of the pulsed component in 2-10 keV en- 
ergy band (Saito 1998). 

B1821-24 From ASCA observations in 0.7-10 keV 

band, its pulsed emission subtracted the background 
and the steady components can be fitted by a power- 
law with a = — 1.2 (Saito 1997). The unabsorbed flux 
in this energy range is 3.5 x 10~"'^^erg cm~^s~"'^, which 
leads to iVpi = 1.4 x 10"d2(ro/roLc)"^. 
B0656+14 Combining ROSAT and ASCA data, 
Greiveldinger ct al. (1996) reported that the X- 
ray spectrum consists of three components: the soft 
surface blackbody with kTg = 67 eV and Ag = 
4.5A*(d/0.76)2, a hard blackbody with kT^ = 129 
eV and Ah = 3.2 x 10-2^^(^/0.76)2, and a power 
law with a = -1.5 and A/pi = 3.1 x lO^cP. The hard 
blackbody component takes the major role in main- 
taining the gap, by virtue of its large emitting area. 
Geminga The X-ray spectrum consists of two com- 
ponents: the soft surface blackbody with kTg — 50 eV 
and As = 0.22A*((i/0.16)2 and a hard power law with 
a = -1.6 and A^pi = 3.9 x 10*^^^ (Halpern & Wang 
1997). A parallax distance of 160pc was estimated 
from HST observations (Caraveo et al. 1996). 
B1055-52 Combining ROSAT and ASCA data, 
Greiveldinger et al. (1996) reported that the X-ray 
spectrum consists of two components: a soft black- 
body with kTs = 68 eV and A^ = 7.3^,(^/1.53)2 
and a hard blackbody with fcTh = 320 eV and A^ = 
2.3 X 10-4^(^/1.53)2. 

J0437-4715 Using ROSAT and EUVE data (Becker 

& Triimper 1993; Halpern ct at. 1996), Zavlin and 
Pavlov (1998) demonstrated that both the spectra 
and the light curves of its soft X-ray radiation can 
originate from hot polar caps with a nonuniform tem- 
perature distribution and be modeled by a step-like 
functions having two different temperatures. The 
first component is the emission from heated polar- 
cap core with temperature kT^ = 10^-^^ K mea- 
sured at the surface and with an area Ah = 5.3 x 
10^'* A* (d/0. 180)2. The second one can be interpreted 
as a cooler rim around the polar cap on the neutron 
star surface with temperature kTg = 10^-^^ K and 
with an area Ag = 0.16A*(d/0. 180)2. Considering 
the gravitational redshift factor of 0.76, the best-fit 
temperatures observed at infinity become kTg = 22 
eV and kT^ = 95 eV. From parallax measurements, 
its distance is reported to be 178 ± 26 pc (Sandhu et 
al. 1997). We adopt d = 180 pc as a representative 
value. 
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4.2. Input Infrared Field 

We next consider the infrared photon field illumi- 
nating the gap. In addition to the references cited 
below, see also Thompson et al. (1999) for Crab, 
B1509-58, Vela, B195H-32, Geminga, and B1055-52. 

Crab Interpolating the phase-averaged color spec- 
trum in UV, U, B, V, R (Percival et al. 1993), J, H, 
K (Eikenberry et al. 1997) bands, and the radio ob- 
servation at 8.4 GHz (MofFctt and Hankins 1996), the 
spectrum in IR energy range can be fitted by a single 
power-law (fig. ^) 

vF^ = 5.75 X lO^^zyi-i^jy (41) 

where v is in Hz. The Crab pulsar's flux at 8.4 GHz 
is, in fact, very uncertain, because it can scintillate 
away from Earth for tens of minutes to hours. Mof- 
fett obtained a value of 0.61 mJy from the profiles he 
collected, while Frail got a value of 0.5 ±0.1 mJy from 
VLA imaging (Moffett 2000, private communication). 
They are within one sigma (0.1 mJy) of each other. 
To estimate a conservative infrared photon number, 
we simply assume that the flux at 8.4 GHz is 0.6 ±0.2 
mJy as a crude estimate, because Moffett (1996) pre- 
viously gave the value of 0.76 mJy. If the error bar at 
8.4 GHz is smaller, then the infrared photon density 
around 0.01 eV reduces further; this in turn results in 
a less upscattered, TeV flux. Equation (^ gives the 
following infrared number spectrum: 

dNm/dei^ = 1.5 x 10"d2(^p/^^^)-2gj^-o.88^ (42) 

where eYRin^c^ refers to the infrared photon energy. 

B0540-69 Its de-extincted optical and soft X- 
ray pulsed flux densities can be interpolated as = 
0.207i^~"-3 Jy (Middleditch & Pennypacker 1985). 
This line extrapolates to 0.47 mJy at 640 MHz, which 
is consistent with an observed value 0.4 mJy (Manch- 
ester et al. 1993). We thus extrapolates the relation 
Fi, = 0.207i^^" '^ to the infrared energies and obtain 

dNi^/dem - 1.6 x lO^^d^ro/mi^cy^ei^-^-'' (43) 

B1509-58 The flux densities emitted from close 
to the neutron star in radio (Taylor, Manchester, & 
Lyne 1993), optical (Caraveo, Mereghetti, & Bignami 
1994), soft X-ray (Seward et al. 1984), and hard X- 
ray (Kawai et al. 1993) bands can be fitted by a single 
power-law F^ = 1.36i^~° '^^ Jy. We thus adopt 

dNi^/dem - 4.7 x 10"d2(ro/n7Lc)"'eiR"'-' (44) 



Vela The flux densities emitted from close to 
the neutron star in radio (Taylor, Manchester, & 
Lyne 1993; Downs, Reichley, & Morris 1973), opti- 
cal (Manchester et al. 1980) bands can be fitted by 
F^ = 1.71 X 10'^z/-°-9i Jy. We thus adopt 

dNm/dti^, = 1.8 X 10^d2(^^/^^^)-2^^^-i.9 (45) 

B1951+32, The flux densities emitted from 
close to the neutron star in radio (Taylor, Manch- 
ester, & Lyne 1993) and soft X-ray (Safi-Harb & 
Ogelman, & Finley 1995) bands can be interpolated 
as Fy = 32.8i^~° '^9 Jy. We thus adopt 

dNi^/dem = 6.3 x 10i°d2(ro/tnLc)"'eiR~' ' (46) 

B0656+14 The nonthermal (most likely, mag- 
netospheric) emission in I, R, V, B bands shows spec- 
tral index —1.53 (Caraveo et al. 1994a; Kurt et 
al. 1998), which is much softer than that of Crab 
(eq.|^). If we were to extrapolate it to the radio 
frequency, the flux density at 1 GHz would exceed 
100 Jy; therefore, it is likely that the soft power-law 
spectrum becomes harder below a cirtain frequency. 
We thus estimate the upper limit of infrared flux den- 
sity by interpolating between radio (6 and 4 mJy at 
0.4 and 1.4 GHz, respectively, Taylor, Manchester, 
& Lyne 1993) and the infrared-optical observations 
(0.71, 0.47, 0.40 iiiy at I, R, V bands). The result is 

dNin/dem = 2.0 x lO^d^ro/wLcyhm-^-^ (47) 

Geminga The upper limit of the flux density in 
radio band (Taylor, Manchester, & Lyne 1993) and 
the flux density in optical band (Shearer et al. 1998) 
gives spectral index greater (or harder) than —0.69. 
Interpolating the infrared flux with these two frequen- 
cies, we obtain 

dNm/dem = 1.9 x I0^d^{r„/wi^cyhm-'-' , (48) 

which gives a conservative upper limit of infrared pho- 
ton number density under the assumption of a single 
power-law interpolation. 

B1055-52 The flux densities emitted from close 
to the neutron star in radio (Taylor, Manchester, & 
Lyne 1993) and optical (Mignani et al. 1997) bands 
can be interpolated as F^, = 4.2 x lO^iy^^ '^'^ Jy. We 
thus adopt 

dNm/dem = 7.7 x Wd^iro/vui^cyhm-'-^ (49) 
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J1617-5055, J0822-4300, B1821-24, and J0437- 
4715 There have been no available infrarad or op- 
tical observations for these four pulsars. We thus sim- 
ply assume that a — —1.5 for these four pulsars and 
that i^F^ = 10^ JyHz at 0.01 eV. We then obtain 

dNin/dem = 2.5 x 10^d2(^^/^^^)-2^j^-i.5 
for J1617-5055, while 

dNm/dem, = 2.1 x 10''d^{ro/wi^c)-^em-'-^ (51) 
for J0822-4300, 

dNm/dem = 1.3 x lO'^'d^ro/wi^cyhiK-^-^ (52) 
for B1821-24, and 

dNm/dem = 3.6 x lO'^d^iro/mi^cy^em^^-' , (53) 
for J0437-4715. 

4.3. Curvature-radiated 7-rays 

In this subsection, we present the results of GeV 
emission via curvature process and compare them 
with observations. The results are presented for the 
two assumed values of a = 30° and 45°. 



For Crab, B0540-69, B1509-58, and J1617-5055, 
the X-ray field is dominated by the power-law compo- 
nent to have high number densities above lO^^cm^'^. 
In this case, the gap half widths are less than 10% of 
tULc- The intrinsic luminosities of these young pul- 
sars in GeV energy range exceed lO'^^ergs s~^. Except 
for the distant pulsar B0540-69, their GeV fluxes are 
expected to be large enough to be observed with a 
space 7-ray telescope. 

For the relatively young pulsars J0822-4300 and 
Vela, the X-ray field is dominated by the surface 
blackbody component; the number density (iVx) be- 
comes about lO^^cm^'^. The gap half width is about 
5% of wlc a-nd the intrinsic GeV luminosity is ^ 
lO^^ergs s^^ 

For the middle-aged pulsar B1951-I-32, its rela- 
tively strong magnetic field at the gap center {B ^ 
lO^G) results in a strong GeV emission like the young 
pulsars. 

The millisecond pulsar B1821-24 has a very strong 
magnetic field {B ^ 10^'^G) at the gap center. How- 
ever, its strong {N^ ~ lO^^'^cm""^) X-ray field pre- 
vents the gap to extend in the magnetosphere. As a 



result, Lqcy is relatively small compared with other 
pulsars. 

For the three middle-aged pulsars B0656-f 14, Geminga, 
and B1055-52, their power-law components are too 
weak to dominate the surface blackbody emissions; 
the surface emissions are also weak (iVx < lO^'''^) to 
allow the gaps to extend more than 10% of njLc- How- 
ever, the extended gaps do not mean large intrinsic 
GeV luminosities, because their small magnetic fields 
{B < lO'^'^G) suppress the acceleration field. In the 
case of Geminga, its proximity leads to a large GeV 
fiux. 

In the case of the millisecond pulsar J0437-4715, 
its weak X-ray field {N^ < lO^^cm"^) due to the soft 
and hard blackbody emissions results in an extended 
gap. This active gap, together with its proximity, 
leads to a large GeV flux next to Geminga. How- 
ever, its relatively strong magnetic field {B ^ 10^ 
G) at the gap center may indicate the presence of an 
additional power-law component, which reduces the 
gap width and hence the GeV flux. Therefore, fur- 
ther hard X-ray observations are necessary for this 
millisecond pulsar. 

4.4. Invisibility of TeV pulses 

If an electron or a positron is migrating with 
Lorcntz factor F ^ 1 in an isotropic photon field, 
it upscatters the soft photons to produce the follow- 
ing number spectrum of 7-rays (Blumcnthal & Gould 
1970): 



c diViR dem 



<PN 3 

— — (Tt 

dtd€j 4 F^ deiR eiR 



2qlTLq+{l + 2q){l-q) + 



{Qq?{l-q) 



2(1 + Qg) 



,(54) 



where Q = 4eiRr and q = €~^/Q{T — e-^); here, e-y 
refers to the energy of the upscattered photons in 
rricC^ unit. Substituting the infrared photon spectrum 
dNiR,/ deiR,, integrating d'^N/dtde-y over eiR, and mul- 
tiplying the 7-ray energy (e^mcC^) and the electron 
number (iVc) in the gap, we obtain the flux density of 
the upscattered, TeV photons as a function of e^. 

We compute the TeV spectra of individual pulsars 
and summarize the results in table 3; the uF^ peak 
frequencies and the fluxes are given in the last two 
columns. Moreover, for the three brightest pulsars 
(Crab, B0656-M4, B1509-58), their computed vF^ 
spectra are presented in flgure ^. It follows from 
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Table 3: Expected 7-ray properties 



pulsar 


a-i 






Ec 




Lgcv/(P 


^pcak 




)TeV* 




deg 






GcV 


erg 




Jy-Hz 


Hz 


JyHz 


Crab 


30 


0.82 


0.023 


4.2 


1.3 • 




2.1 


• 10^^ 


2.5 


• 1027 


1.9 • 


IOI2 




45 


0.96 


0.016 


5.2 


1.8 • 


1033 


3.0 


• 10^2 


2.7 


• 102'^ 


4.4 • 


10^2 


B0540-69 


30 


0.21 


0.048 


7.8 


8.5 • 


1033 


3.6 


•lOio 


2.7 


•1027 


8.3 • 


lOi" 




45 


0.24 


0.034 


9.8 


1.3- 


1034 


5.3 


•101° 


2.9 


•1027 


2.1 • 


10" 


B1509-58 


30 


0.21 


0.075 


5.1 


5.7- 


1033 


3.1 


•IO12 


3.4 


•1027 


1.1 • 


10" 




45 


0.23 


0.053 


6.6 


9.0- 


1033 


4.9 


•1012 


3.6 


•1027 


2.9 • 


10" 


J1617-5055 


30 


0.093 


0.097 


8.7 


1.4- 


1034 


7.4 


•1012 


1.2 


•1027 


1.0 


108 




45 


0.10 


0.068 


11.1 


2.2- 


10^4 


1.1 


•10^3 


2.0 


•1027 


2.5 


108 


J0822-4300 


30 


0.69 


0.045 


2.5 


5.7- 


1032 


1.2 


•1012 


6.0 


• 1026 


1.0 


107 




45 


0.57 


0.035 


3.9 


1.4- 


1033 


3.1 


•1012 


7.6 


.1026 


3.7 


107 


Vela 


30 


0.27 


0.084 


3.5 


1.7- 


1033 


7.0 


•10^3 


3.2 


.1026 


8.4 • 


lO^o 




45 


0.26 


0.062 


4.9 


3.3 • 


1033 


1.4 


•1014 


4.0 


• 1026 


2.6 • 


10" 


B1951+32 


30 


0.16 


0.092 


5.4 


2.8- 


1033 


4.7 


•1012 


6.3 


.1026 


3.1 • 


101° 




45 


0.18 


0.064 


6.9 


4.2- 


1033 


7.1 


•10^2 


7.4 


.1026 


8.0 • 


10^° 


B1821-24 


30 


0.22 


0.067 


5.2 


4.6 • 


1032 


1.8 


•10" 


5.0 


• 1026 


6.2 


107 




45 


0.25 


0.047 


6.7 


7.2- 


1032 


2.9 


•10" 


4.7 


• 1026 


1.6 


108 


B0656+14 


30 


0.43 


0.160 


1.1 


2.0- 


1032 


3.7 


•10^2 


3.7 


.1026 


1.7 • 


10^2 




45 


0.69 


0.098 


1.2 


1.8- 


1032 


3.2 


•1012 


3.8 


• 1026 


3.0 • 


1012 


Geminga 


30 


0.059 


0.334 


4.0 


4.1 • 


1033 


1.7 


•1015 


1.5 


•1027 


7.0 


109 




45 


0.083 


0.213 


4.5 


4.3- 


1033 


1.8 


•10^5 


6.5 


.1026 


1.4 • 


10^° 


B1055-52 


30 


0.62 


0.130 


0.98 


8.5 • 


10^1 


3.8 


•10" 


6.0 


• 1026 


6.3 


109 




45 


i.ot 


0.078 


1.0 


7.1 • 


1031 


3.2 


•10" 


6.0 


•1026 


1.1 • 


101° 


J0437-4715 


30 


0.043 


0.371 


5.0 


6.1 • 


1032 


2.0 


•10" 


6.3 


.1026 


5.1 


107 




45 


0.062 


0.233 


5.4 


6.1 • 


1032 


2.0 


•10" 


5.9 


.1026 


9.9 


107 



^ The Lorentz factors are assumed to be the value that can be attained if particles are 

accelerated by in length H. 

t The TeV flux are evcilucLtecl at the lyFi/ peak frequency, i^peak* 
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table 3 and figure || that the TeV emissions are in- 
visible with the current ground-based telescopes, ex- 
cept for Crab and B0656-I-14. Since H becomes as 
small as Zacc for the Crab pulsar, significant fraction 
of the particles are, in fact, unsaturated. Therefore, 
the expected TeV flux is overestimated. To constrain 
the absolute pulsed TeV flux from the Crab pulsar, 
we must discard the mono-energetic approximation 
for the particle distribution function and explicitly 
solve the Boltzmann equations of particles and 7- 
rays, together with the Poisson equation for <I>, un- 
der suitable boundary conditions. In addition, con- 
sidering the fact that the tertiary infrared photons 
are not isotropic but have small collision angles with 
the particles, we can understand that the TeV flux 
computed from equation (JsJ) are, in general, overes- 
timated. For B1055-52 and B0656-hl4, the emitting 
areas of the soft blackbody components are unnat- 
urally large. Therefore, more accurate X-ray obser- 
vations are necessary for a quantitative prediction of 
their TeV fluxes. 

It is worth noting that the TeV fluxes are less than 
3% of the GeV fluxes except for Crab, B0540-69, and 
B0656-I-14; this conclusion is qualitatively consistent 
with Romani (1996). The predicted TeV flux from 
B0540-69 is greater than the GeV flux, because (for 
ai — 30° for instance) its infrared photon number 
density in the energy interval 0.01 eV and 0.1 eV at- 
tains 7.5 X lO^^cm"'^, which well exceeds the X-ray 
number density 6.4 x IQ^^cm"'^ between 0.1 keV and 
100 keV. 

5. Discussion 
5.1. Summary 

To sum up, we have considered the electrodynamic 
structure of an outer gap accelerator in which rel- 
ativistic particles emit 7-rays via curvature process. 
Imposing the gap closure condition that a single pair 
produces one pair in the gap on average, we solve 
selfconsistently the gap width as a function of the X- 
ray flelds and the pulsar parameters. Once the gap 
width is known, we can further compute the acceler- 
ation fleld and the resultant 7-ray emissions. It was 
demonstrated that the luminosities of GeV and TeV 
emissions are a decreasing function of the X-ray en- 
ergy and number density. We also showed that the 
expected vF^ fluxes (< 10^^'^Jy Hz) of IC-scattered, 
TeV 7-rays from the outer gaps of rotation-powered 
pulsars are less than the observational upper limits. 




frequency [Hz] 



Fig. 5. — A single power-law fit of phase-averaged 
color spectrum of the Crab pulsar. The abscissa is 
the photon frequency in Hz, while the ordinate is the 
photon flux in Jy-Hz. 




10"" 10" 
frequency [Hz] 



Fig. 6. — Expected TeV spectra of the three bright- 
est pulsars: Crab (solid lines), B0656-I-14 (dashed 
ones), and B1509-58 (dotted ones). The thick and 
thin curves represent the cases of ai = 30° and 45°, 
respectively. The abscissa is the photon frequency in 
Hz, while the ordinate is the photon flux in Jy-Hz. 
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except for Crab and B0656+14. For Crab, energy- 
dependent particle distribution function should be 
considered, whereas for B0656+14, more accurate X- 
ray observations are required. It is concluded that the 
difficulty of excessive TeV emission, which appears in 
the CHR picture, does not arise in the present outer 
gap model. 

5.2. Stability of the Gap 

The outer gap in the present model is stable, re- 
gardless of whether the X-ray field is dominated by 
a surface blackbody or a magnetospheric power-law 
component. Consider the case when the gap width 
H slightly increases as an initial perturbation. It in- 
creases both E\\ and V^ap, which in turn increases 
both Ec and iV^. The increase of results in the 
decrease of Eth- 

(1) When the surface blackbody dominates the X- 
ray field, the X-ray spectrum and luminosity are un- 
changed by the perturbation. Therefore, the decrease 
of Et\i implies the decrease of Ai or A2 and hence Ap. 

(2) When the magnetospheric emission dominates, 
the secondary and tertiary emissions will increase 
with Ec and A^^; therefore, A^pi increases as well. Ac- 
cordingly, the decrease of i?th and the increase of A^pi 
imply a significant decrease of A3 and hence Ap. 

In either case, it follows that Ap decreases owing to 
the initial increase of H . Reminding the gap closure 
condition H = Xp/Nj, we find a negative feedback 
which cancels the initial perturbation of H . 

5.3. Pulse Sharpness 

Let us discuss the expected sharpness of GeV 
pulses. It seems unlikely that the azimuthal width 
of the gap increases with decreasing H. Therefore, 
it would be possible to argue that the solid angle in 
which the primary 7-rays are emitted decreases as the 
arc of the gap along the last open field line (i.e., 2H) 
decreases. On these grounds, we can expect a sharp 
pulse when h 1 holds, such as for Crab and J0822- 
4300. 

Qualitatively speaking, the same conclusion can 
be expected for millisecond pulsars and magnetars. 
In the case of a millisecond pulsar, its fast rotation 
shrinks the light cylinder. In such a small- volume 
magnetosphere, the outer gap is immersed in a dense 
magnetospheric, power-law X-ray emission. As a re- 
sult, Ap decreases to reduce h. In the case of a mag- 
netar, its strong magnetic field makes the expansion 



coefficient A in equation (||) be large. Therefore, a 
very thin (h <^ 1) gap with a strong would be ex- 
pected. In short, for young pulsars, millisecond pul- 
sars, and magnetars, their high-energy pulsations are 
expected to show sharp peaks. 

5.4. Validity of Assumptions 

First, we reduced the Poisson equation into the 
one-dimensional form (eq.[^), by assuming D± ^ H. 
Let us briefly consider the two-dimensional effect due 
to the transfield derivative in the Poisson equation. 
When D± becomes small, the gap shifts outwards, 
£^ll is partially screened, and H enlarges (fig. 12 in 
Paper I; see also Cheng, Ho, & Rudermann 1986a for 
a screened, or spatially constant E'n in a thin gap). 
Owing to the screened acceleration field, the GeV and 
TeV fluxes becomes small compared with those ob- 
tained in D± ^ H case. On these grounds, we can 
constrain the upper limit of the TeV fluxes in the 
transversely thick limit, D± ^ iJ." 

Secondly, let us discuss the case when the assump- 
tion of the vacuum gap breaks down. In this case, 
the charges in the gap partially cancel the original E\\ 
obtained in the vacuum gap (eq. §). The partially 
screened E^^ results in the decrease of the TeV fluxes. 
On these grounds, we can regard the TeV fluxes pre- 
sented in the present paper as the firm upper limits. 

Thirdly, we consider the influence of cyclotron res- 
onance scatterings. For one thing, the soft blackbody 
emission from the whole surface may be scattered to 
be anisotropic (Daugherty & Harding 1989). Such 
effects are important for polar cap models, because 
the collision angles (cos^^ /ic) suffer significant cor- 
rections. Nevertheless, in an outer gap, such correc- 
tions are negligibly small. Moreover, the cyclotron 
resonance increases the effective emitting area and 
decreases the temperature. For simplicity, we neglect 
these two effects in this paper, because they cancel 
each other. For example, the decreased tempera- 
ture results in a decrease of the target photons above 
a certain threshold energy for pair production. On 
the other hand, the increased emitting area increases 
the number of target photons above the threshold, 
thereby cancel the effect of the decreased tempera- 
ture. What is more, the /larrf blackbody emission from 
the heated polar caps may be scattered to be smeared 
out. That is, most of the hard X-rays may be scat- 
tered back to the stellar surface owing to cyclotron 
resonance scatterings and reemitted as soft X-rays 
(Halpern & Ruderman 1993). In this case, the hard 
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component will be indistinguishable with the origi- 
nal soft component due to the neutron-star cooling. 
Nevertheless, for older pulsars such as B0656-I-14 and 
B1055-52, these effect seems to be ineffective proba- 
bly due to their less dense electrons around the polar 
cap near the neutron star surface. 



5.5. Gamma-ray Luminosity vs. 
Luminosity 



Spin-dov^rn 



Curvature-radiated luminosity, iccVi has a weak 
dependence on the spin-down luminosity, igpin, if we 
fix the transfield thickness of the gap, D^/wi^c. In 
another word, the evolution of D^/wi^q is crucial to 
discuss the icov oc ispin°'^ relation (Thompton et al. 
1994; Nel et al. 1996). To solve D±, we must analyze 
the two-dimensional Poisson equation on the poloidal 
plane; however, it is out of the scope of the present 
paper. 

5.6. Synchrotron Radiation Below^ 10 MeV 
Energies 

As we have seen, the accelerated particles reach 
curvature-radiation reaction limit to become roughly 
monoenergetic. The curvature spectrum in lower en- 
ergies then becomes a power law with a spectral in- 
dex 1/3, which is much harder than the observed 7- 
ray pulsar spectra. In this subsection, we demon- 
strate that the 7-ray spectrum below a certain energy 
(say 10 MeV) is dominated by a synchrotron radia- 
tion from freshly born particles and that the expected 
7-ray spectra further softens. 

As an example exhibiting a soft power-law 7-ray 
spectrum from eV to GeV energies, we consider the 
Crab pulsar. To discriminate whether curvature or 
synchrotron process dominates, we separately con- 
sider each process and take the ratio of the radiation- 
reaction forces. That is, we ignore much complicated 
synchro-curvature process, because such details are 
not important for the present purpose. 

Let us first consider the case of — 45°, which 
gives i?5 = 6.5 X 10^ at the gap center. Since the 
curvature-radiated 7-ray energy is = 5.2 GeV, 
the freshly born particles have the Lorentz factors of 
Fq '--^ 5 X 10^. A particle with this Lorentz factor emit 
synchrotron radiation around the energy 



^^sync — 



S/iFo^eSsin Xp 

47rTOcC 



(55) 



We can solve the evolution of the Lorentz factor 
and Xp simultaneously by the method described in 
§ 5.3 of Paper I. We present the evolution of sinxp due 
to synchrotron radiation shortly after the pair produc- 
tion in figure ^ and the evolution of the longitudinal 
momenta in figure ^j. In both figures, the abscissa 
designates the distance along the fieldlines in njLc 
unit with respect to the birth place (distance^ 0). 
The particles are supposed to be created with posi- 
tive momenta; therefore, electrons turn back to have 
negative longitudinal momenta. Positrons lose lon- 
gitudinal momenta on the initial stage of accelera- 
tion, because the relativistic beaming effect causes 
the synchrotron-radiation-reaction force not only in 
the transverse but also in the longitudinal directions. 

It follows from figure that we can approximate 
sinxo ~ 0.3 when the particles have not run A/ ~ 
1 X 10~^n7Lc for ai = 45°. When Xp is kept around 
0.3, Lorentz factors are also in the same order of Fq. 
Substituting sinXp = 0.3 and Fq ~ 5 x lO'^ into equa- 



tion (55), we obtain hi^sync ~ 8.5 MeV as the cen- 
tral energy of the synchrotron spectrum for ai = 45°. 
The fraction of the particles having sinXp ^ 0.3 and 
Lorentz factor ^ Fq to the saturated particles then 
becomes 

AN. Al ^^^^ 



H 



6 X 10" 



We can estimate the ratio between the synchrotron 
radiation from the freshly born particles and the cur- 
vature radiation from the saturated particles as fol- 
lows: 

AN,, dPsync/di^ 



No dPcurv/dl^' 



where 



dP. 



sync 



^/Se^B sinxo 



dv 



F 



rricC 



hv 



sync 



dPc. 



Rc \ En 



dv 



F{x) = X K5/3{y)dy] 



(57) 



(58) 



(59) 



(60) 



7^5/3 refers to the modified Bessel function of 5/3 or- 
der; F in equation ( |5^ ) denotes the saturated Lorentz 
factor and becomes 2.4 x 10'' for 45° for Crab. At the 
synchrotron peak energy, 0.29/ii'sync — 2.5 MeV, the 
ratio becomes i?sc = 12. 

In the same manner we can consider the case of 



where Xp denotes the pitch angle of the particles. 



a-i = 30°. In this case, we have hi>s 



1.7 MeV and 
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ANc/Nc - 1 X 10"''. As a result, we obtain Rgc = 11 
at 0.5 MeV. 

We can therefore conlude that the 7-ray spectrum 
befow certain energy (~ 10 MeV) is dominated by the 
synchrotron radiation from freshly born particles. 

In addition, in the case of Crab, the unsatu- 
rated motion of particles (Zacc ^ H) implies that the 
synchro-curvature radiation from unsaturated parti- 
cles are also important. Therefore, the spectrum 
below GeV will become much softer compared with 
the simple curvature spectrum with central energy 

^ 5.2 (or 4.2) GeV for a-, = 45° (or 30°). 

5.7. Comparison with Zhang and Cheng model 

Finally, we point out the difference between the 
present work and Zhang and Cheng (1997); they con- 
sidered a gap closure condition so that the curvature- 
radiated 7-ray energy may be adjusted just above the 
threshold of pair production. That is, they considered 
the 7-ray energy to be about E^^zc = (TOcC^)^/i?x, 
where Ey^ refers to the characteristic X-ray energy. By 
equating E^^zc with the central energy of curvature 
radiation (eq. jl^ in our notation), they closed the 
equations. When the soft (or hard) blackbody emis- 
sion dominates, Ey^ can be approximated by SfcTg (or 
3/cTh). 

The model of Zhang and Cheng (1997) is, in fact, 
qualitatively consistent with our gap closure condi- 
tion, provided that the X-ray are supplied by the soft 
or hard blackbody emission. More specifically, our 
model gives about 2 times larger characteristic 7-ray 
energy compared with their model. To see this, we 
present in figu re p| the ratio between E,, computed 
from equation (|14D and E^^zc'i the hard blackbody or 
the power-law components are not considered in this 
calculation. The abscissa indicates the soft blackbody 
temperature, fcTg. For the three thick curves, O2 is 
fixed at 0.5; the solid, dashed, and dotted hues cor- 
responds to 1130 = 1.0, 3.0, and 0.3, respectively. For 
the two thin curves, on the other hand, /Z30 is fixed 
at 1.0; the dashed and dotted curves corresponds to 
= 1.0 and 0.25, respectively. At small fcTg, our 
model gives more than twice greater 7-ray energy 
compared with Zhang and Cheng (1997); neverthe- 
less, the difference is not very prominent. 

It should be noted, however, that the spectra of 
the X-ray radiation are explicitly considered in our 
present model in the sense we perform the integra- 
tion over X-ray energies in equations (p^), (pT|), and 




distance from birth point 



Fig. 7. — Pitch angle evolution due to synchrotron 
radiation as a function of the position. The abscissa 
denotes the distance (in tolc unit) along the field lines 
from the birth place (distance= 0). The thick and 
thin lines correspond to ai — 30° and = 45°, re- 
spectively. 
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( P5| ) and that the additional, power-law component is 
considered in our present model. 

It would be interesting to investigate the back re- 
action of the accelerated particles on the polar cap 
heating, which was deeply investigated by Zhang and 
Cheng (1997). Consider the case when the voltage 
drop in the gap approaches the surface EMF, V^,. 
Such an active gap will supply copious relativistic pri- 
mary particles to heat up the polar cap due to bom- 
bardment. The resultant hard blackbody emission 
supplies target photons for pair production to make 
the realistic solution deviate from the thick solid curve 
and approach thin curves at small kTg (fig. ^ . There- 
fore, this sort of back reaction on the X-ray field due 
to the relativistic particles needs further considera- 
tion. 



-a 

3 



2000 



-2000 





y 

w A-- 
\ \ 







-10*^ 10"^ 2*10*^ 

distance from birth point 



Fig. 8. — Longitudinal momentum evolution due to 
synchrotron radiation as a function of the position. 
The ordinate are nomalized in meC^ unit. The ab- 
scissa and the lines are the same as figure 7. 
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Fig. 9. — The ratio between Ec computed from equa- 
tion (13) and -£7,70 (see text). Neither the hard 
blackbody nor the power-law components are consid- 
ered. X-rays are supplied by the whole surface black- 
body emission; the abscissa refers to the temperature, 
kT,. 
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